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Study of Gluon versus Quark Fragmentation in T ^ gg'j and 
e'^e~ qq^ Events at i/s=10 GeV 


Abstract 


Using data collected with the CLEO II detector at the Cornell Electron 
Storage Ring, we determine the ratio -Rchrg for the mean charged multiplicity 
observed in T(1S) —> gg^ events, < >, to the mean charged multiplicity 

observed in e+e“ ^ qq'j events, < >. We find Rchrg = = 

^^quark^ 

1.04 ± 0.02(stat) ± 0.05(syst) for jet-jet masses less than 7 GeV. 
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I. INTRODUCTION 


Understanding hadronization, or the process by which elementary partons (gluons and 
quarks) evolve into mesons and baryons, is complicated by its intrinsically non-perturbative 
nature. Naively, one expects that because of the greater color charge of gluons compared 
to quarks, radiation of secondary and tertiary gluons is more likely when hadronization is 
initiated by a gluon rather than a quark. This results in a greater number of hnal state 
hadrons as well as a larger average opening angle between the hadrons in the former case 
compared to the latter case. In the limit > cxo, the ratio of the number of hadrons 
produced in gluon-initiated jets to the number of hadrons produced in quark-initiated jets 
is expected, in lowest order, to approach the color degeneracy factor of 9/4 [Q]. 

Many experiments have searched for, and found, multiplicity and jet shape differences 
between quark and gluon fragmentation [0]- p^. At energies, qqg events are readily 


distinguished by their three-jet topology. Within such events, quark and gluon jets can be 
separated by a variety of techniques including vertex tagging. Because gluons rarely fragment 
into heavy quarks, they will produce jets that form a vertex at the e+e" interaction point. 
Quark jets, to the contrary, tend to form a detached vertex when the jet contains a long-lived 
bottom or charm quark. Unfortunately, the assignment of final state hadrons to the initial 
state partons is rarely unambiguous and relies on Monte Carlo simulations to determine the 
fraction of times that an observed hadron is correctly traced to a primary parton. 

The 10 GeV center of mass energy range offers a unique opportunity to probe quark and 
gluon fragmentation effects, without relying on Monte Carlo simulation to associate the final 
state hadrons with an initial state parton. The decay T(IS') — gg'y allows one to compare 
the gg system in a gg^y event with the qq system in e’''e“ —»■ qq'f events. In both cases, the 
system recoiling against the photon consists (to lowest order) of hadrons that have evolved 
from either a two-gluon or a quark-antiquark system. The properties of the recoil systems 
can then be compared directly.^ 


II. DETECTOR AND DATA SAMPLE 


The CLEO II detector is a general purpose solenoidal magnet spectrometer and 
calorimeter. The detector was designed for efficient triggering and reconstruction of two- 
photon, tan-pair, and hadronic events. Measurements of charged particle momenta are made 
with three nested coaxial drift chambers consisting of 6, 10, and 51 layers, respectively. These 
chambers fill the volume from r=3 cm to r=l m, with r the radial coordinate relative to the 
beam (z) axis. This system is very efficient (e >98%) for detecting tracks that have transverse 
momenta {pt) relative to the beam axis greater than 200 MeV/c, and that are contained 
within the good fiducial volume of the drift chamber {\cos6\ <0.94, with 6 defined as the 


^Although there may be gluon radiation from the initial partons, we do not distinguish such 
radiation explicitly in this analysis. Thus, the states that we are comparing are, strictly speaking, 
gg'^ and qqj to lowest-order only; additional gluon radiation, to which we are not experimentally 
sensitive, may be present in many of the events in our sample. 
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polar angle relative to the beam axis). Below this threshold, the charged particle detection 
efficiency in the fidncial volnme decreases to approximately 90% at pt ~100 MeV/c. For 
Pt <100 MeV/c, the efficiency decreases ronghly linearly to zero at a threshold of pr ~30 
MeV/c. 

Beyond the time of flight system is the electromagnetic calorimeter, consisting of 7800 
thallinm doped Csl crystals. The central “barrel” region of the calorimeter covers abont 
75% of the solid angle and has an energy resolntion of 


E 


(%) 


0.35 


+ 1.9-0.IE; 


( 1 ) 


E is the shower energy in GeV. This parameterization translates to an energy resolntion of 
abont 4% at 100 MeV and 1.2% at 5 GeV. Two end-cap regions of the crystal calorimeter 
extend solid angle coverage to abont 95% of dvr, althongh energy resolntion is not as good as 
that of the barrel region. The tracking system, time of flight connters, and calorimeter are 
all contained within a 1.5 Tesla snpercondncting coil. Flnx retnrn and tracking chambers 
nsed for mnon detection are located immediately ontside the coil and in the two end-cap 
regions. 

We nse 63 pb“^ of data collected at the T(1S) resonance (y^=9.46 GeV) as a sonrce of gg'y 
events and 198 pb“^ at y/s=10.52 GeV (on the continnnm jnst below the T(4S) resonance) as 
a sonrce of qq'j events. The 7 in onr qq'j sample resnlts primarily from initial state radiation 
(ISR) [^]. We compare events for which the invariant masses of the gg and qq systems 
recoiling against the hard photon {Mrecou, defined by Mrecou = - E^/Ebeam)) are 

the same. 

To snppress low-mnltiplicity QED events, we reqnire that the thrnst of the event (cal- 
cnlated nsing all the photon candidates and good qnality charged tracks) be less than 0.97, 
and that there be at least 3 good charged tracks in the event. Photon candidates are se¬ 
lected from showers with widths and patterns of energy deposition consistent with that of 
a photon, as opposed to nentral hadrons (e.g. merged vr^’s, E)), nentrons, etc.). Owing to 
the excellent ability of the GLEO-II detector to distingnish high-energy tt^’s from photons, 
approximately 75% of the potential 7 r° backgronnd is removed at this stage by shower topol¬ 
ogy cnts. Additional 7 r° snppression is achieved with an explicit 7 r° —>• 77 mass cnt, to be 
discussed later. To ensure that the events are well-contained within the GLEO detector, we 
require | cos 6 '..y| < 0.75 ( 6 *.y is defined as before as the polar angle between the beam axis and 
the direct photon). 

In addition to the gg'j and ggq samples, we search for T(1S)—>• ggg and e+e“ ^ qq events 
containing a tt^ of energy comparable to the photon in the gamma-tagged sample. These 
samples, referred to as ggiii^) and q'g(vr°) are used to quantify the background levels from 
7 r° decays. 


III. EXPERIMENTAL TECHNIQUE 

One of the most basic parameters used to characterize any event is the mean charged 
track multiplicity, < Nchrg >■ We plot < Nchrg > as a function of the mass recoiling against 
the direct photon, covering the recoil mass range from 4 to 7 GeV. This mass interval 
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corresponds to 0.45 < E^/E^ea.ra < 0.82 for the gg'-f sample. We require the recoil mass to be 
greater than 4 GeV to ensure that we are signihcantly above the ccq threshold; we expect 
that our gg jet sample therefore includes u,d,s,c quarks approximately uniformly, apart from 
small phase space effects.^ 

To determine the characteristics of ggy events, we must subtract the background from 
non-resonant ggy and > rry events produced in annihilations at a/s = Mxis- 

This can be done by direct scaling of the event sample collected at the continuum of the 
T(4S) resonance. In addition to ggy events, this latter sample of events also includes rry 
events, as well as any other continuum backgrounds which may be present in our resonant 
sample. By subtracting this “raw” sample from our T(1S) sample, we therefore account 
for all of the non-resonant backgrounds to ggy at the T(1S) energy. We hnd that before 
subtraction, this background to ggy comprises about 10% of the gamma-tagged sample taken 
at the T(1S) energy. 

Similarly, in order to isolate ggy events at 10.52 GeV, we must quantify rry contami¬ 
nation. This is done using a Monte Garlo simulation of tan pair events. We hnd that rry 
events comprise about 10% of the ggy data sample passing the other event selection cuts 
specihed above. Beam gas and two photon backgrounds were investigated and found to be 
negligibly small. 

To determine the level of gg(vr°) contamination to ggj and qq{7i^) contamination to ggy, 
we hrst determine how often a single high energy photon can be matched with other photons 
in the same event to form a 7r°. Such photons are explicitly vetoed as likely tt^ daughters. 
The 7r° contamination, of course, results from cases for which we hnd only one of the two 
7 r° daughter photons in the detector [^. From Monte Garlo simulations, we determine the 
probability that the second 7r° daughter photon will be found. The likelihood of detecting the 
second photon varies from 50% at Mrecou = 6.5 GeV to 77% at Myecou = 4 GeV. Knowing 
the fraction of times that the second daughter photon goes undetected, and the total number 
of TT^’s that we reconstruct in our sample, we thereby determine the fraction of our direct 
photon candidates which are actually 7r° daughters, but are not vetoed as such. This is 
shown in (Figure 2) as a function of recoil mass. 

Our backgrounds are summarized below: 


Sample 

Vi 

Background (%) 

991 

9.46 GeV 

Sum of all non-T(lS) events (10%) 

991 


T(1S)—gg(7 r°) contamination to T(1S)—>• gg'j (~ 5%) 

qqi 

10.52 GeV 

rry (10%) 

qqi 


—>• gg(7r°) contamination to ggy (~ 8%) 

qqi 


77 (< 1%) 

qqi 


beam-gas (< 1%) 


^We note that the qq mixture in our qqy sample may be different for cases in which the photon is 
emitted in the final state as opposed to the initial state due to the different diagrams responsible 
for these processes. We expect our sample to be dominated by ISR events, as is evidenced by the 
photon polar angle distribution. 
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FIG. 1. vr^ contamination, as a function of recoil mass, in gg^ events (squares), and in qq^ 
events (circles); the y axis gives the fraction of our apparent gg^ (or qq^) event sample which are 
actually gg{'K^) (or qq{'K^)) events. 
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To determine how our measured charged multiplicity values are biased by the 7r° back¬ 
ground, we make a comparison plot of the charged multiplicity distribution for gg'y vs. 
gg{7^^) events, and a similar plot for ggy vs. qq{T^^) events. From Figure |^, we see that the 
charged multiplicity distributions for the gamma-tagged and the vr'^-tagged samples are sim- 
ilar.0 Quantitatively, the effect of 7r° backgrounds is to reduce the measured value of Rchrg 
by about 1% relative to the true value. We statistically correct for this effect in subsequent 
plots, and include the uncertainty due to 7r° backgrounds in our overall systematic error. 

In Figure we show the multiplicity distributions in bins of recoil mass for our 
background-corrected ggj and qq'j samples. From the distributions in this Figure, we extract 
the mean multiplicities as a function of gg or qq mass, for pure gg'j and qq'y samples. As 
shown in Figure the ratio of < Nchrg > resulting from gluon fragmentation to < N^hrg > 
from quark fragmentation is Rchrg = 1-04 ± 0.02, after all the aforementioned background 
corrections. 


IV. SYSTEMATIC ERRORS 

In order to investigate the sensitivity of the result to the Nchrg > 3 requirement and 
to suppress rry backgrounds, we compare the ratios obtained from Nchrg > 3 with ratios 
obtained using the tightened multiplicity interval 5< Nchrg <11- This comparison also 
partially addresses possible biases in our measured ratio due to acceptance effects as well 
as the effect of the minimum charged multiplicity cut. If the shapes of the true charged 
multiplicity distributions were much different for quarks vs. gluons at low multiplicities (for 
which our acceptance is worse than at high multiplicity, due to the effect of our hadronic 
event selection cuts as well as trigger effects), then our derived average will be artihcially 
‘pulled’ by our Nchrg >3 requirement. A complete Monte Carlo simulation study of possible 
differences between the detector level multiplicities and the true multiplicities, including 
unfolding effects, indicates negligible bias. The value of the ratio using Nchrg >5 is found to be 
statistically consistent with that for Nchrg >3. The effect of imposition of the minimum Nchrg 
cut has also been checked by repeating the entire analysis, but using the neutral multiplicity 
Nneut rather than the charged multiplicity as the comparison variable. For this cross-check, 
we retain the cut Nchrg > 3, but make no cut on the minimum neutral multiplicity. Before 
any consideration of backgrounds, this cross-check yields Rneut = 1.06T0.03 (statistical error 
only). This value is in good agreement with our measured value of Rchrg- 

Although we compare systems of equivalent mass, the energies of the systems recoiling 
against the photon are higher for qq events than for gg events at an equivalent invariant 
mass. On average, the ggy system will therefore be more collimated. We investigate possible 
energy dependent effects with Monte Carlo simulation, by looking at systems of equal recoil 
mass, but generated at different energies (A/i=10.52 GeV vs. A/i=9.46 GeV). We hnd a 


^It is of some interest to note the similarity in multiplicity between the gg'^ and the gg{TT^) 
samples. Although the first process contains two primary gluons, whereas the second contains 
three primary gluons, this plot qualitatively suggests that at these energy scales the available 
energy for fragmentation is the main factor in determining the final state charged multiplicity. 







nts 


2400 


3031296-014 


1600 


4 GeV < M < 5 GeV 
recoil 


997 


997 




Charged Multiplicity Ratio for 997 /qq 7 


1.10 


3031296-015 



FIG. 4. Charged multiplicity ratio for background subtracted gg^/qq^ data. 
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systematic shift of 0.006 in the ratios. We can also evaluate energy effects by direct inspection 
of Figure |. Although the Rchrg values measured for the lowest recoil mass are statistically 
consistent with the Rchrg values measured for the highest recoil mass considered, our results 
are also consistent with Rchrg increasing slightly with recoil mass. We conservatively include 
a systematic error of 2 % to account for any possible systematic variation. 

We also compare the gg'^/qq'y multiplicity ratios from data taken during three distinct 
running periods, covering different running conditions and triggers. The separate data sets 
give ratios which are all statistically consistent with each other. Nevertheless, we conser¬ 
vatively include the maximum variation of 2 % observed between the three data sets in the 
overall systematic error. 

Finally, we consider the dependence of the mean multiplicity on the dip angle, | cos 6 '..),|, of 
the direct photon. We expect that initial state radiation qq'y events tend to produce photons 
that are more forward peaked (have larger values of | cos0.^|) than for gg'j. For qq'y events, 
we therefore expect that the particles emerging on the other side of the event, opposite the 
high energy 7 , are more likely to be lost down the beam pipe. We have therefore explicitly 
considered any possible dependence of Rchrg on \cos9^\ by plotting Rchrg in different bins of 
I cos 6 *..^ I . At our level of statistical precision, we do not observe any such systematic effect 
and conclude that the measured Rchrg value is not measurably biased by such acceptance 
effects. 0 We explicitly quantify this by measuring i?chrg in two regions of cos6^\ |cos 6 '.y| <0.35, 
and 0.35 < \cos9^\ <0.70. The variation between the two ratios is found to be 0.02 ± 0.03. 
We conservatively assign the error on this value as a measure of the maximum possible 
systematic error due to such polar angle acceptance effects. 

In summary, the systematic errors are: 


Systematic Error 

Value 

7 r° background 

0.003 

Recoil mass-dependent effects 

0.006 

Truncated (5< Nchrg <11) vs. non-truncated mean 

0.02 

Variation between data sets 

0.02 

Rchrg dependence on Mj-econ 

0.02 

cos 6 * cut dependence 

0.03 

Total 

0.05 


V. SUMMARY 

The ratio of <Nchrg> for gluons to <Nchrg> for quarks measured here is smaller than 
those found by the OPAL, ALEPH, SLD and DELPHI experiments, at ^/s ~ M^o. The 


“^There are potentially two competing effects here - although the photon is more forward peaked 
for qq'y events, it is also the case that quark jets should be more narrowly collimated than gg^ 
events. This may tend to mitigate any effects due to the difference in photon angular distributions, 
since charged hadrons produced in qq'y events should be more contained within the fiducial region 
of our drift chamber. 
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ratios compare as follows: 


Gollaboration 

<N>g / <N>g 

Kinematic Regime 

GLEO 96 

1.04 ±0.05 

< Ejet >< 7 GeV 

DELPHI 

11 

1.24 ±0.015 

< Ejet > = 10-40 GeV 

DELPHI 

11 

1.06 ±0.18 

< Ejet > = 10 GeV 

SLD p 


1.36 ±0.24 

< Ejet > = 24 GeV 

OPAL 93 

18| 

1.27 ±0.04 

< Ejet > = 24 GeV 

ALEPH [| 


1.19 ±0.04 

< Ejet > = 24 GeV 

OPAL 96 

21 * 

1.58 ±0.03 

< Ejet > = 39 GeV 


*Tliis OPAL result is for < N >g —1/ < N >g —1, considering only uds quarks 


Fuster and Marti |^, DELPHI [Q] and Fodor |Q have recently considered how the ratio 
of charged multiplicity for gluon to quark jets depends on the hadronic center-of-mass en¬ 
ergy. Phenomenologically, as the energy scale increases, the likelihood of successive partons 
radiating also increases. Having a greater color degeneracy than quarks, the multiplicity of 
gluon-initiated jets increases faster than the multiplicity of quark-initiated jets. This causes 
Rchrg to increase slowly with energy. Our result is smaller than the LEP results for Rchrg, 


showing the expected energy dependence 
surement by DELPHI at ^/s ~10 GeV, albeit with substantially better precision 


Our result is also consistent with the mea- 

We note 


that the qq mixture in this experiment is likely to be different than the quark mixture for the 
high energy experiments mentioned above. Recent analyses by OPAL and ALEPH 


in fact, have explicitly measured the difference in ratio values when comparing gluon jets to 
b-quark jets versus gluon jets to uds-quark jets. 
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